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ABSTRACT 

Formulas  and  tables  are  developed  for  compressing 
the  power  series  representation  of  a  function  to  obtain 
minimal  order  approximating  polynomials  maintaining 
specified  accuracy  throughout  an  assigned  range  of 
arguments.  Coefficients  for  polynomials  approximat¬ 
ing  the  function  1-e^  are  derived  and  scaled  for  use  in 
a  fixed  point  computer. 
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I. 


INTRODUCTION 


In  certain  launch  vehicle  guidance  equations  (Ref.  1)  there  occurs  a 
function  of  the  form 


u(v)  =  1  -  e'^  (1) 

Values  of  the  argument  lie  within  a  range  -0.  3  <  v  <  2.  1.  It  is  required  that 

any  approximation  used  uhall  yield  u  =  0  when  v  =  0.  The  accuracy  desired 
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is  such  that  the  leaat  significant  bit  of  u  should  be  2 

At  present,  u(v)  is  approximated  by  an  eighth  order  polynomial  derived 
by  truncating  the  Maclaurin  series  for  the  exponential. 


i=l 


(2) 


The  approximation  (2),  however,  provides  the  desired  accuracy  only  when 

V  <  1.  At  the  upper  end  of  the  range  of  arguments,  the  significance  deteri- 
-8 

orates  to  about  2 

This  report  gives  the  derivation  of  coefficients  for  two  new  approximat¬ 
ing  polynomials.  One  of  these  is,  and  the  other  maybe,  sufficiently  accurate 
over  the  entire  pertinent  range  of  arguments.  Section  II  describes  a  general 
method  by  which  a  long  power  series  can  be  compressed  to  yield  the  approxi¬ 
mating  polynomial  of  minimum  order  to  maintain  specified  accuracy  throughout 
a  given  range  of  arguments.  In  Section  III,  the  method  is  applied  to  the  function 
defined  in  Equation  (1).  Section  IV  introduces  considerations  of  scaling,  neces¬ 
sary  to  retain  significance  while  avoiding  overflows  in  a  fixed  point  computer, 
and  convenient  to  enable  bit-by-bit  simulation  of  the  computer  operations  for 
purposes  of  checking  and  error  analysis.  In  Section  V,  coefficients  are  evalu¬ 
ated  for  both  seventh  and  eighth  order  polynomials  to  approximate  function  (1) 
over  a  limited  range  of  arguments.  The  errors  of  these  approximations  are 
compared  in  Section  VI  with  those  of  the  currently  used  approximation  (2).  In 
Section  VII  it  is  shown  that  the  approximating  polynomials  derived  in  this  report 
are  preferable  to  Equation  (2)  for  numerical  differentiation  of  Equation  (1). 


II.  COMPRESSION  OF  POWER  SERIES 


In  Reference  2  it  is  shown  that  although  a  function  may  be  expressed  by  a 
convergent  infinite  power  series,  there  are  other  expressions  more  suitable  for 
approximating  the  function  within  specified  tolerance  for  a  limited  range  of  the 
argument.  Of  particular  importance  are  the  Chebyshev  polynomials  T.(x)  and 
Tii'(x).  The  former  are  useful  when  -1  <  x  <  1 ,  and  the  latter  when  0  <  x  <  1. 
Both  have  the  property  that,  within  the  range  of  applicability,  their  maxima 
are  +1  auid  their  minima  are  -I.  A  series  of  Chebyshev  polynomials  converges 
more  rapidly  than  any  other  series  representing  a  function  when  arguments  are 
restricted  to  the  appropriate  range. 

In  this  report,  the  polynomials  T^(x)  are  used,  with  0  <  x<,  1.  The 
normalization  of  the  original  independent  variable  to  this  range  is  discussed 
in  Section  III. 


T.(X)  = 

j=0 

where 

t  _  1 

‘oo  -  7 

tjo  U>1] 

'ij 


(4) 


In  Reference  2,  the  definition  of  Tq(x)  is  given  special  attention.  For  present 
purposes,  however,  the  value  of  1/2  is  appropriate.  Table  II-l  lists  t^j  for 
0  <  j  <  i  <  7. 


It  it  further  shown  in  Reference  2  that  any  integral  power  of  x  can  be 
expressed  as  a  finite  aeries  of  Chebythev  polynomials. 


where 


'(nl”)  t0< 


Table  II-2  lists  hj,(n)  for  0  i  <  13  and  7  ^n<  13. 


If  m  is  chosen  <  n,  then  substituting  (3)  into  (5)  for  0  <  i  <  m  -  1  gives 


n  ^  2-(2n-l)  te  gj(m,n)  x^  +  ^h.(n)T» 


T?(x) 


where 


g.(m,n)  =  t.  .  [0  <  j  <  m  -  1  <  n] 

i=j 


gj(7,n)  and  gj(8,n)  are  listed  in  Tables  II-3  and  II-4  for  n  <  13. 

Equation  (7)  shows  how  any  power  n  ^  m  of  x  can  be  represented  by  a 
polynomial  in  x  of  order  m  -  1  and  a  linear  combination  of  T?(x)  with 
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m  <  i  <  n.  The  coefficients  gj(m,n) 
observed  that 


and  h^(n)  are  all  integers. 


It  should  be 


gj(m,n)  + 


^\.(n)  =  2^"-^ 

i=m 


If  0  <  X  <  1.  then  the  series 


(9) 


can  be  expressed  as  the  sum  of  a  polynomial  in  x,  of  some  order  m  -  1,  and 
a  residual  series  of  Chebyshev  polynomials  of  order  >  m. 


Sm-1  QQ 

c.x*=  >^d  .X*  +  ® 

1=0  i=m 


:Tf{x) 


(10) 


From  (7), 


00  m- 1 


oo 


Zc  x'  -  Y^c  x‘  + 

i  /  .1  i  /  ^  n 


i=0  i=0 


n=m 


’m- 1 


n 


^^g.{m,  n)x^  +  y^i(n)Tf(x) 


i=0 


i=m 


(11) 


Equating  coefficients  of  and  T:^(x)  in  (10)  and  (11)  gives 


oo 


^  =  c.  +  ^^2"^^”"^^c^g.(m,n)  [0<i<m-l] 


mi 


n=m 


(12) 


e. 

1 


(13) 


=  [m  <  i  <  n] 

n=i 

Equation  (13)  provides  a  basis  for  selecting  the  order,  m  -  1,  of  an  approxi¬ 
mating  polynomial  whose  error  shall  not  exceed  a  specified  magnitude  when 
0  <  X  <  1.  Since,  in  this  range,  |T?(x)  |  £  the  error  of  approximation  is 
limited  by 


i=m  n=m  i=m 

By  trial  and  error,  a  minimal  m  can  be  found  such  that  E^,  computed  by 
(14),  is  sufficiently  small.  With  m  thus  determined,  the  coefficients  of  the 
approximating  polynomial  cam  be  found  from  (12). 

In  practice,  the  infinite  sum  in  (12)  can  be  replaced  by  a  finite  sum.  If 
the  power  series 


r 

converges  when  0  <  x<  1,  then  for  some  r  the  terms  beyond  c^x  can  be 
neglected.  Equation  (10)  then  becomes 


m-1 


oo 


i=0  i=m  i=r+l 


(15) 
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The  error  of  approximation  it  limited,  for  0  <  x  <  1 ,  by 


(16) 


The  m  coefficients  of  the  approximating  polynomial  are  then 


d'  . 
mi 


n) 


n=m 


[0<i<m  -  1] 


(17) 


-7- 


-8- 


i 


{ 

i 

I 


33,649  -3,083,472  45,224,256  -241,196,032 


177,100  -16,096,100  234,416,000  -1,237,716,480 


f 
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III.  APPLICATION  TO  u(v)  =  1  -  e'' 

In  order  to  apply  the  method  developed  above  to  the  function  defined  in 
equation  (1),  it  is  first  necessary  to  transform  the  variable  to  one  lying  in  the 
range  between  0  and  1.  If  b  <  v  <  a  +  b,  then 


Since  it  is  required  that  the  polynomial  approximation  u  (v)  vanish  when 
V  =  0,  and  since  it  is  further  desirable  that  the  approximation  be  most  accurate 
in  the  neighborhood  of  v  s  0,  the  process  of  Section  II  is  applied  to  the  function 


uW  .  1  -  ^  b)"  . 


V  -(ax  +  b)  7  j  {n  +  1)1 

n-0  i=( 


'c.x*  [0  <  x<  1] 

\  rnmm 


Expanding  (ax  +  b)”  by  the  binomial  theorem,  it  follows  that 


.  oo 


-  a  bJ 

i  “  TT^^i  +  j  +  l)(j! 
j  =  0 


A  value  of  r  is  easily  chosen  such  that 


oo 


Ek' 


i=r+l 


(19) 


(20) 


can  be  neglected  for  b  ^  v  ^  a  +  b.  Application  of  equation  (16)  leads  to  a 

selection  of  m  which  keeps  vE  within  tolerance,  m  is  the  number  of 

mr 


-13- 


coefficients  in  the  approximating  polynomial.  Then,  using  (17)  to  determine 


d'  .. 
mi 


m-1 


m-1 


u  (v) 
m' 


-V 


'k  .vJ 
mj 


(21) 


i=0 


From 


i 

(v  -  b)^ 

j=0 

It  follows  that 

m-1 

i=j 

Note  that  if  b  =  0,  then  c.  =  a'/(i  +  1)1  and  k  .  =  a”'^d'  .. 

X  '  mj  mj 

It  must  be  emphasized  that  the  accuracy  of  approximating  u(v)  by 

deteriorates  very  rapidly  when  v  exceeds  the  range  whose  limits  have  been 

used  in  determining  the  coefficients  k  .. 

*  mj 


mi 


[0  <  j  <  m  -  1] 


(22) 


IV.  COMPUTER  SCALING 


The  polynomials  (v)  are  to  be  evaluated  in  a  fractional  binary 
computer  with  word  length  of  23  magnitude  bits  and  a  sign  bit.  A  rounded 
multiply  operation  is  available. 

In  programming  this  computer,  it  is  customary  to  assign  a  scale 
factor  s(N)  to  a  number  N  such  that  |n|  <  2*.  In  a  word  representing  N, 
the  value  of  the  least  significant  bit  is  2*'^^.  Numbers  to  be  added  must  be 
scaled  alike.  The  scale  of  a  product  is  the  sum  of  the  scales  of  the  factors. 

In  order  to  realize  the  advantage  of  rounded  multiplication,  there 
should  be  no  shifting  of  intermediate  products  in  a  program  for  polynomial 
evaluation.  It  follows  from  (21)  and  the  arithmetic  of  scale  factors  that 


Maximum  precision  and  the  avoidance  of  overflow  in  intermediate 
computations  require,  in  the  present  application,  that 


s(v)  =  2 
s(u^)  =  5 

8(k^.)  =  3  -  2j 


(24) 


The  bit  configuration  in  a  computer  word  for  a  number  N  with  a  scale 

8  8  I 

factor  8  is  the  same  as  that  for  the  number  2  N  with  a  scale  factor  s  . 
Instead  of  computing  k^^  from  (22)  and  assigning  individual  scale  factors 
from  (24),  it  is  preferable  to  transform  the  polynomial  coefficients  to  integers 
which  can  all  be  scaled  with  s  =  23.  The  advantage  lies  in  the  greater 
convenience  with  which  the  computer  operations,  including  round-off,  can  be 
simulated  bit  by  bit  for  purposes  of  error  analysis  and  program  checking. 


Let 


Then 


where 


K  .  =  2^°‘'‘^jk  .. 
mj  mj 


TT 

U  =2  u  , 
m  m 


^  “  fc» 

W  =  2  V, 


s(K  .)  =  23  1 
'  mj'  ' 


8(U  )  =  23 
'  m' 


s(w)  =  0 


(25) 


m-  i 


U  (w)  =  -  w  ^  K  . 
m'  mj 


j=0 


(26) 


U_  with  s  =  23  is  the  same  configuration  of  digits  as  u__  with  s  =  5. 
m  "  ^  m 
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V.  EVALUATION  OF  COEFFICIENTS 

The  range  -0. 3  <  v  <  2. 1  includes  values  beyond  all  those  to  be 
expected  in  the  current  application.  Therefore,  a  =  2.4;  b  =  -0.  3.  With 
these  parameters,  the  c^  are  computed  by  equation  (20)  and  listed  in 
Table  V-1  for  0  <1  <  17. 

«  1  fi 

Accuracy  to  the  order  of  2“  is  desired  in  the  approximation  to  u(v) . 

Since 

2.  1^  Cj  <  2"^^  , 

i=14 


while 


GO 


2. 1^^  c.  is  nearly  2 
1=13 
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the  value  of  r  =  13  is  used  for  termination  of  the  sums  in  equations  (16)  and 
(17). 

From  (16), 


2. 1  E.^  j3  =  754-  10’®  <  2“^^ 

2.  1  Eg  =  78-10"®  <  2"^° 

It  is  evident  that  m  =  8  keeps  the  error  well  below  the  specified  tolerance, 
while  m  =  7  may  be  acceptable  if  slightly  larger  errors  can  be  tolerated 
when  V  >  1.  The  difference  between  r  +  1  =  14  and  m  =  7  or  8  measures 
the  benefit  gained  by  the  series  compression  method  of  Section  II  in 
minimizing  the  order  of  an  approximating  polynomial  satisfactory  throughout 
the  range  of  arguments. 
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/ 


Values  of  computed  by  (17),  with  m  =  7  and  m  s  8,  are  listed 
in  Table  V>2. 

From  (25)  and  (22),  with  a  =  2.4,  b  =  -0.3, 


m-1 

=  2^°"^^^  (0.3)"-^  [0<j<m-l] 

i=j 


(27) 


Values  of  computed  by  (27),  with  m  s  7  and  m  =  8  are  listed  in 
Table  V-3. 


Table  V-1. 


c.  from  Equation  (20)  with  a  =  2.  4, 


b  =  -0.3 


i 

c. 

1 

0 

.86393  92644 

1 

.98496  83497 

2 

.76789  18791 

3 

.45365  10976 

4 

.21551  84202 

5 

.08557  59883 

6 

.02917  93574 

7 

.00871  64991 

8 

.00231  64848 

i 

c. 

1 

9 

.00055  44093 

10 

.00012  06818 

11 

.00002  40891 

12 

.00000  44398 

13 

.00000  07600 

14 

.00000  01214 

15 

.00000  00182 

16 

.00000  00026 

17 

.00000  00003 

% 
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Table  V-2.  from  Equation  (17)  with  r  =  13,  m  =  7  and  8 


i 

•^Si 

<^71 

0 

.86393  90612 

. 86394 

22776 

1 

.98499  39577 

. 98467 

87486 

2 

.76736  44347 

.77240 

77891 

3 

.45776  39967 

.42750 

38699 

4 

.19992  35730 

.28638 

10778 

5 

.11713  94957 

.00966 

48447 

6 

-.00501  19813 

.08720 

93573 

7 

.02634  89539 

Table  V-3.  K _ .  from  Equation  (27)  with  m  =  7  zmd  8 

rnj 


j 

*^8j 

K7. 

0 

1,048,576 

3,048,577 

1 

2,097, 148 

2.097,259 

2 

2,796, 363 

2,795,331 

3 

2,796,485 

2,779.653 

4 

2,  214,973 

2,433,572 

5 

1,645,497 

751,671 

6 

405,  518 

1.959,997 

7 

986,971 
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VI.  ERROR  ANALYSIS 


L> 


Calculations  of  the  seventh  order  polynomial  U^,  equation  (Z6),  were 

made  with  from  Table  V-3.  All  multiplications  were  rounded,  to  simulate 

computer  operation.  The  differences  between  these  results  and  2^®  (l-e'^)  are 

plotted  in  Figure  VI-1.  When  -0.  36  v  <  1,  the  errors  lie  virtually  within  the 

shaded  region  of  round-off  error,  ±  1/2  •  2"^®.  When  1  ^  v  ^  2.  1,  the 

- 1 8 

oscillatory  error  function  never  exceeds  2.  5  ■  2  in  magnitude. 

In  contrast,  the  eighth  order  truncated  Maclaurin  series  approximation 
departs  rapidly  from  u{v)  when  v  >  1.  The  magnitude  of  its  error  reaches 
2"^®  when  v  =  1.2,  2”^^  when  v  =  1.5,  and  2  when  v  =  1.9. 

Errors  of  the  eighth  order  approximating  polynomial  Ug  are  wholly 
obscured  by  round-off,  and  therefore  are  not  shown  in  Figure  VI- 1. 

The  warning  is  repeated  that  approximating  polynomials  derived  by 
the  methods  of  this  report  cannot  be  expected  to  remain  valid  outside  the 
range  of  arguments  for  which  they  are  specified. 
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VII.  NUMERICAL  DIFFERENTIATION 


Reference  1  recommends  that,  since  du/dt  is  used  in  the  application, 
u  should  be  approximated  by  truncation  of  the  exponential  series.  However, 
the  polynomials  derived  in  this  report  improve  the  approximation  to  the 
derivative  as  well  as  the  approximation  to  the  function  itself. 

If  an  approximation  to  u(v)  is  denoted  by  u’^(v)  and  the  error  by  c  (v),  then 

u(v^)  =  u^lv.)  +  t(v^)  (28) 

u(v2)  -  u(Vj)  =  u’f'{v2)  -  u^lvj)  +  -  €(Vj)]  (29) 

In  words,  if  the  differentiation  is  to  be  approximated  by  differencing, 
then  the  error  of  a  difference  is  given  by  the  difference  of  the  errors. 

It  is  evident  from  Figure  VI- 1  that  in  the  region  1  <  v  <  2.  1,  the  slope 
of  the  error  curve  for  the  truncated  Maclaurin  series  is  significantly  greater 
in  magnitude  than  the  slope  of  the  error  curve  for  u.y(v).  Therefore,  even 
the  seventh  order  polynomial  derived  here  gives  a  better  approximation  for 
the  derivative  than  does  the  series  truncated  to  eighth  order. 
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